Shear-induced overaging in a polymer glass 
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A phenomenon recently coined as "overaging" implies a slowdown in the collective (slow) relax- 
ation modes of a glass when a transient shear strain is imposed. We are able to reproduce this 
behavior in simulations of a supercooled polymer melt by imposing instantaneous shear deforma- 
tions. The increases in relaxation times Ari/2 rise rapidly with deformation, becoming exponential 
in the plastic regime. This "overaging" is distinct from standard aging. We find increases in pres- 
sure, bond-orientational order and in the average energy of the inherent structures (< eis >) of the 
system, all dependent on the size of the deformation. The observed change in behavior from elastic 
to plastic deformation suggests a link to the physics of the "jammed state". 

PACS numbers: 61.43.Fs, 61.25.Hq, 62.20.Fe, 64.70.Pf 



In recent years, there has been considerable progress 
geared towards understanding how glasses respond to 
shear. Phenomena such as shear thinning and "rejuve- 
nation" are common when shear flow is imposed. Un- 
like crystals, glasses "aga," meaning that their state de- 
pends on their history 0, |^ • When a glass falls out 
of equilibrium, it evolves over very long timescales. It 
has been found that relaxation in supercooled liquids of- 
ten depends on cooperative and spatially correlated mo- 
tion, meaning the dynamics are heterogeneous (J| . Much 
research has also been devoted to developing a wide- 
reachingjCoherent theory which can explain the jammed 
state jj, Q . It has been found that such a state is charac- 
terized by the appearance of a yield shear stress and can 
be achieved by changing the load, density or temperature 
of a system. In addition, it has been found that the con- 
cept of random close packing is somewhat ill-defined and 
in the so-called jammed state, one can increase the degree 
packing at the expense of randomness, or vice- versa, thus 
allowing for a variety of possible jammed configurations 

Ii,8j. 

The concept of jammed state also provides insight as 
to the effect of stress on a glass. When imposing a con- 
tinuous shear strain, one can rejuvenate the glass, ef- 
fectively wiping out the memory of the system llOjl . 
In this sense, the glassy state is particularly sensitive to 
shearing. In 2002, Viasnoff and Lequeux published ex- 
perimental results for colloidal suspensions, showing that 
after imposing transient shear, one obtains simultaneous 
overaging and rejuvenation, since the relaxation times are 
altered in a non-trivial manner In this case, over- 

aging means that relaxation times of the system become 
longer more quickly than is normally the case. Their 
results are qualitatively explained through Bouchaud's 
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Trap Model and the related Soft Glassy Rheology (SGR) 
model [lol IT^ . In essence, this phenomenological ap- 
proach relies on the distribution of relaxation times r 
(or equivalently, of potential wells of different depths) 
which can be over- or underpopulated when applying a 
transient shear or a temperature step. Using a similar 
perspective. Lacks and Osborne have argued, in a zero 
temperature study, that the so-called "overaging" is dif- 
ferent in nature from ordinary aging, in that the minima 
visited can be similar in both cases, but not identical 
It has also become clear that it is difficult to induce de- 
creases in energy by strain in well-annealed glasses, since 
the disappearance of energy minima does not usually lead 
the system to a lower potential well 0, l2l| . 

In this letter, we report on some of the characteris- 
tics of this overaging from a more general perspective. 
We find that overaging is present and easily observed in 
simulations of a common glass-forming polymer model 
simply by imposing relatively small, instantaneous shear 
deformations. After a certain waiting time t^, one can 
observe an unambiguous slowing-down in the decorrela- 
tion of particle positions, just as was seen experimentally 
We use the term "overaging" to refer only to the 
longer relaxation times, not to any decrease in energy 
of the system (as in the case of ordinary aging). We 
have identified two distinct regimes of overaging, cor- 
responding to elastic and plastic deformations. There 
is a rapid increase in relaxation times Ari/2(e), which 
contrasts with the associated increase in the average en- 
ergy of the internal structure of the system (< e/s >), 
as introduced by Stillinger and Weber We concur 

with the recent study, finding that overaging is quite dis- 
tinct from aging in its regular sense [1^]. In our case, 
we show that it is possible to obtain longer relaxation 
times without the system's slow evolution towards more 
energetically favorable configurations. In addition, the 
yield shear strain of the material plays a key role in this 
phenomenon. Finally, we resolve an increase in order in 



2 



the system, primarily associated with the elastic energy 
stored in the system. 

The simulations are performed by molecular dynam- 
ics, using the velocity Verlet algorithm, combined with 
velocity-rescaling to achieve constant temperature po|. 
We adopt the bead-spring model to simulate the poly- 
mer melt 'l6], each chain consisting of 10 monomers 
linked together through a finitely-extensible nonlinear 
elastic (FENE) potential of the form UpENEifij) — 

0.5fci?§log l-{rij/Raf , where k = SOeij/cr^ and 

Rq = 1.5(7. In addition, all particles interact via 
a truncated and shifted Lennard- Jones (LJ) potential, 

with the cutoff 



(cr/ry) 



ULj{nj) = 4£lj 

radius set at 2.5a. For simplicity, all data are presented 
in reduced LJ units based on the mass m, the radius of 
each particle a and the LJ energy scale elj- Each MD 
step corresponds to 0.005 reduced time units. The sam- 
ples have 105 chains, so a total of N = 1050 particles, and 
periodic boundary conditions (PBCs) are applied in all 
directions. We use up to 40 independent samples in order 
to average the results. The combination of the FENE 
and LJ potentials causes two competing length scales, 
thus inhibiting crystallization and producing a typical 
fragile glass JJ, 18, 19] . In a previous paper, we have 
identified a glass transition (GT) based on the simula- 
tion timescales used, as well as a rigidity transition (RT) 
located just below the GT We also observed that 

the shear modulus /x, computed via instantaneous de- 
formations, is highly dependent on the size of the shear 
deformation e. Specifically, larger values of e can be char- 
acterized as irreversible or plastic deformations which al- 
ter the structure (and energy landscape) of the system. 
In such cases, the residual stress is reduced leading to 
smaller fi values. In contrast to the recent work of Lacks 
and Osborne Q|, we have a well-annealed system, as 
well as a non-zero temperature. This means that the 
system has the ability to explore a range of energy wells 
(as opposed to a single one) through thermal activation. 

All samples are set at T = 0.44 and are initially at iden- 
tical pressures in the neighborhood of 0.57 in the rigid 
pha se ( near the onset of rigidity which occurs below the 
GT [17|). They are created via a slow compression of the 
simulation box, followed by a damped-force algorithm to 
realize the NPT ensemble |23|- By definition, these sys- 
tems are out of equilibrium, but we start with samples 
that have been allowed to evolve considerably and see 
no evidence of aging on the timescales of our computer 
experiments. Starting with a cube of side L, an affine 
shear deformation Sxy is applied in the x direction, in a 
plane with its normal along the y direction. An atom ini- 
tially in position (x, y, z) is displaced to (x + txyy, y, z). 
The boundaries of the simulation box consequently are 
shifted for Xmin from to exyU , and for Xmax from L to 
^xyU + L. For a given system, the shearing is repeated 
in the five other directions, substituting xy by yz,—xy, 
etc. Individual samples lack symmetry and therefore the 
various deformations will not usually give the same stress 
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FIG. 1: Correlation function of Eq. showing the effects 
of repeated shear, and different waiting times after each one, 
on the slow relaxation modes of the system, compared with 
the reference (unsheared) case (bold line). Legend for both 
the inset and main frame: first deformation &t = Q (solid 
line) and = 10^ (dashed); second deformation at t^, = 
(dotted) and = 10^ (dash-dot). Main frame: large shear 
deformations (e = 0.2). Note the initial combination of re- 
juvenation and overaging, followed by pure overaging. Inset: 
Small shear (e = 0.05). The relaxation changes, although in 
a very different manner, and causes small overaging. 



components. Once the system has been allowed to relax 
for a fairly long time , this process can be repeated on 
the previously deformed sample, either deforming it fur- 
ther or returning it to its original shape. Both methods 
yield identical increases in relaxation times. 

In order to monitor the relaxation in the system, we 
compute the two-time, "self" part of the incoherent scat- 
tering function. 



Cq(ijj 



1 ^ 
N ^ 



exp[iq-(rj(t^, +T)-rj(i^))] (1) 



where the wavevector q is close to the first peak in the 
structure factor S'(q), is the time elapsed since the 
deformation. This "two-time" correlation function has 
been proven useful in investigating aging, because of the 
different timescales in relaxation 2, 3J. We use ti/2, the 
time it takes for Cq to decrease to 0.5, in order to gauge 
the local decorrelation in the system. The effect of the 
deformations is seen in Fig. ^ As was found in Ref. 
|Tll |. for large deformations (e > 0.1), there is an initial 
combination of overaging and rejuvenation, as relaxation 
times corresponding to both high and low-energy states 
are overpopulated. Eventually, the initial shape of the 
function is recovered since only the slow relaxations re- 
main overpopulated. Interestingly, one can repeat this 
process a number of times, achieving the same results af- 
ter every applied deformation. Our focus is not on the 
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FIG. 2: Energy of the inherent structure (IS) of the sys- 
tem, calculated for various e immediately after (squares) and 
tnj = 10^ (circles) after the deformation. Note that, for large 
(plastic) deformations, the inherent structure of the system is 
able to evolve. 



transient effects, but rather on the behavior after a rel- 
atively long tw As seen in the inset of Fig. ^ small 
shear also causes overaging but is the result of a simpler 
physical mechanism. 

We start from a well-relaxed system. The strain in- 
creases the energy and appears to raise the energy bar- 
riers between available energy wells, as evidenced by the 
increased decorrelation times. For small shear as seen 
in Fig. 121 there is no major configurational change. The 
difference between the < e/5 > immediately after and a 
relatively long time after a given deformation is minimal. 
This difference increases rapidly as we enter the plas- 
tic regime where a large part of the strain energy is re- 
laxed by the configurational changes. Plastic shear maps 
the system onto a very different energy well, from which 
the system can escape fairly easily. This is the source 
of the initial combination of overaging and rejuvenation 
observed in Fig. ^and in Ref. 11]. After a reasonable 
waiting time (f^ = 10'^), there is a net increase in ri/2- 
The characteristic incremental relaxation time Ati 12(e) 
is shown for various shear deformations e in Fig. y] af- 
ter subtracting the reference relaxation time Ti/2(e = 0), 
which is a constant on our simulation timescales. In the 
plastic regime, there is a clear exponential behavior. 

To get some inkling of what may be causing the over- 
aging we look at the structural changes occurring. When 
applying a relatively large deformation, the chains are 
initially stretched from their equilibrium lengths. But 
after letting the system relax for approximately 10'^ time 
units, the radial distribution function and the average 
radius of gyration of the chains return to their "unde- 
formed" values. A more sensitive measure of structural 
change is the local bond-orientational order parameter 
Qe. local which uses higher-order spherical harmonics and 



FIG. 3: Incremental relaxation times after tw ~ W for vari- 
ous e. The line is a guide to the eye and shows an exponential 
dependence on e for plastic deformations. 
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FIG. 4: Evolution of the local bond-orientational order pa- 
rameter Qe, local at tw = 10"^ after the deformation. While the 
statistics are poor, there is a clear increase in Qe, local, at least 
for e < 0.1. Data based on 40 samples deformed twice, and 
treating each deformation as independent. 
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Indeed, Qe, local increases slightly as we increase e, while 
e < 0.1, followed by a plateau or slight decrease (see Fig. 
0J|, noting that subsequent deformations do not have a 
large effect on the order parameter. This ordering, is 
most apparent for small shear, indicating that it is di- 
rectly related to the elastic energy stored in the system. 
Increased order is usually associated with a lower energy 
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FIG. 5: Pressure in the system, calculated for various e at 
tu, — 10^ after the deformation. The pressure only increases 
for small (reversible) deformations. 

state. In our polymeric system, however, the potential 
energy is dominated by the stiff intra-chain bonds. The 
increased order is at the expense of the potential energy, 
and leads to a pressure increase (see Fig. [SJ (also ob- 
served in LJ binary mixtures [2lj|). The increased Ti/2 
is an indication that there are larger energy barriers be- 
tween configurations, i.e. that the strain "jams" the sys- 
tem; and in the case of plastic shear, reorganizes the 
potential wells. The dynamics are also affected. Once 
again, the two regimes of shear strain are very distinct. 
Small strains have little or no effect, while large strains 
reduce the "heterogeneity" in the system (i.e.; there are 
less deviations from a Gaussian distribution of particle 
displacements and more particles become "mobile" 



|24|. This, in turn, causes a decrease in the shear modulus 
/i of the system, i. e. shear softening |l7l | . Increased mo- 
bility is not incompatible with increased relaxation time. 
Average waiting times are dominated by the contribu- 
tions from the largest Ti/2 of the "metabasins" whereas 
the mobility is controlled only by the most mobile parti- 
cles mug. 

An interesting aspect of this study is the similarity 
with "tapping" experiments on granular matter and 
glassy systems \27^. which can produce jammed struc- 
tures. Common features in the structure and force dis- 
tributions of granular an d "g lassy" jamming have also 
recently become apparent |23|. Just as tapping granular 
matter can both increase or decrease the density, tran- 
sient or instantaneous shear can produce a variety of ef- 
fects on the glassy system, while continuous shear always 
"unjams" or rejuvenates it. In conclusion, this is not ag- 
ing in the conventional sense (the system has not aged 
more quickly); < ejs > does not decrease and neither 
does the mobility The very similar behavior observed 
in colloids 0, , the present polymer glass 0| , and 
binary Lennard Jones mixtures 2l| suggests a univer- 
sal behavior and unique origin of the overaging, observed 
so far in two of these systems. Shear produces changes 
in the microstructure which lead to increased relaxation 
times. This may be put to profitable use in some appli- 
cations. It appears also to be a way to increase ordering. 
Finally, repeated applications of the deformations yield 
increasingly larger relaxation times, and may generate a 
glassy state of practical interest, whose properties are yet 
to be fully understood. 
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